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Abstract. In muonic atoms the Uehling potential (an effect of a free electronic vacuum polarization loop) 
is responsible for the leading contribution to the Lamb shift causing the splitting of states with An — 
and Al 7^ 0. Here we consider the Lamb shift in the leading nonrelativistic approximation, i.e., within 
an approach based on a certain Schrodinger equation. That is valid for low and medium Z as long as 
(Za) 2 <C 1. The result is a function of a few parameters, including k — Zam,j,/m e , n and I. We present 
various asymptotics and in particular we study a region of validity of asymptotics with large and small k. 
Special attention is paid to circular states, which are considered in a limit of n> 1. 



PACS. 36.10.Gv Mesonic atoms and molecules, hyperonic atoms and molecules 
and quantum electrodynamic effects in atoms and molecules 



1 Introduction 

The gross structure of energy levels in all kinds of hydrogen- 
like atoms is generally of the same form determined by the 
Schrodinger-Coulomb equation 
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where m is the mass of the orbiting particle which is an 
electron in a conventional atom and a heavier particle in a 
muonic or exotic atom. However, details of the spectrum 
and, in particular, the structure of the energy levels with 
the same value of the principal quantum number n are dif- 
ferent in different kinds of atoms. For example, in muonic 
atoms at low and medium Z the largest splitting between 
states with the same n is the one for states with Al ^ 
(the Lamb splitting) which is essentially a nonrelativistic 
effect. 

In the nonrelativistic approximation the leading con- 
tribution to the Lamb shift in muonic atoms (i. e., the 
Uehling correction) has been known analytically for a while 
for certain levels £Q, however, only numerical results used 
to be quoted in the textbooks (see, e.g., 0). A reason for 
that is the complicated form of the analytic expressions. 
For instance, in the simplest case of the ground state the 
result is of the form J] 
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and m e is the electron mass. Expressions for other states 
are similar, but more complicated. They involve functions 
A(n n ) with a characteristic parameter 
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and coefficients similar to those in Eq. QJ depend on values 
of the principal and orbital quantum numbers, n and I. 

The mass of the orbiting particle m in a non-conventional 
hydrogen-like atom is much above the electron mass m e . 
We consider here the vacuum polarization effects for a 
hydrogen-like atoms with an orbiting particle, which in 
particular may be a muon (m^ ~ 207 m e ; n ~ 1.5 Z), a 
pion (ttt-tt ~ 273 m e ; k ~ 2 Z), an antiproton (m-p ~ 1836 m e ; 
K ~ 13 Z) etc. The relativistic effects for those atoms are 
quite different for various reasons, while the result in the 
leading nonrelativistic approximation is the same. Further 
we do not distinguish between various possibilities of the 
orbiting particles and mainly speak about a muon, but the 
equations could be applied to any orbiting particle. 
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Analytic results have been known for some time even 
for hydrogen-like atoms with a Dirac particle [Hj and since 
recently for the case of a Klein-Gordon particle 0|. They 
are rather cumbersome, containing the hypergeometric func- 
tion 3F2 and far from being transparent. For instance, the 
relativistic result [30 for the nl states reads as a finite 
sum over basic integrals 



n dependence and a determination of a real parameter of 
expansion are important to find the range of validity of 
various asymptotics. 

Here we derive a general expression for the Lamb shift 
at medium values of the nuclear charge Z. Finally the 
vacuum polarization correction is presented in the leading 
nonrelativistic approximation in the form 



K abc (K n ) = -K c n B(a + 1/2, 1 - 6/2 + c/2) 

x 3 F 2 (c/2, c/2 + 1/2, 1 - 6/2 + c/2; 

1/2, a + 3/2- 6/2 + c/2; K n ) (5) 

- °- k^ +1 B(a + 1/2, 3/2 - 6/2 + c/2) 

x 3^2(0/2 + 1, c/2 + 1/2, 3/2 - 6/2 + c/2; 
3/2, a+ 2 -6/2 + c/2; n 2 n ) , 

where ai^ftti, ck2, «3; /?i , /?2 ; z) stands for the generalized 
hypergeometric function (see, e.g., and B(ai, ot2) is the 
beta function. The parameters a, 6, c are linear functions 
of n and I in the nonrelativistic case, while for the rela- 
tivistic results they contain certain additions of relativistic 
corrections which go to zero at the limit of (Za) — > 0. The 
argument of 3F2, K n , is reduced to K 2 n in the nonrelativistic 
approximation. 

Meanwhile, muonic atoms offer a special region of pa- 
rameters where the result can be essentially simplified (see 
e.g. 013)- For instance, the Uehling correction for the 
ground state [S||3] (cf. Eq.QJ) takes the form 
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The Lamb shift splits the levels with An = and Al ^ 
and for this reason we also consider a specific difference 



AE(ls) ~ -- (Zafmc 2 

IT 



2 , x 11 



(6) 



@nii'{K n ) = F n i(n n ) - F nV (n n ) 



(9) 



The simplification is possible because in the range of medium 
Z we can apply for the ground state a double expansion 
over two parameters: 



and typically for our calculations V — I + 1. 

We find in this paper asymptotics for low and high n n 
and determine regions of their validity. We study in more 
detail circular states and show that for them the 1ow-k„ 
expansion is an expansion over n ■ K n = k, while the high- 
n n asymptotics is actually an expansion over n 2 / k. 

Additionally to well-defined regions of these expan- 
sions (k < 1 or n 2 J ' k <C 1) there are also two intermediate 
regions: 

— low K n , when K n ^C 1, but n ■ K n ~ 1; 

— high K n , when K n ~^> 1, but K n /n ~ 1. 

We discuss behavior of the Uehling correction in these two 
specific regions. 



2 The Uehling correction in the 
nonrelativistic approximation: general 
consideration for the Lamb shift 

Let us first remind how the Uehling correction is calcu- 
lated in a general case. The Lamb shift in muonic atoms 
is a result of perturbing the Coulomb potential 



Za< 1 , 

K> 1 . 
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Here and further we consider only a leading non-relativstic 
approximation (i.e., the leading term of the Za expan- 
sion). 

Highly excited states in muonic and exotic atoms are 
of particular interest for precision measurements because 
they offer a certain suppression of the interaction between 
the nucleus and the orbiting particle. The n dependence of 
theoretical expressions, even of the simplest asymptotics, 
is not a trivial issue. One can see from expressions with 
the generalized hypergeometric function 3.F2 that while 
the argument is «„, the parameters are n dependent and 
in fact in actual situations some are proportional to n. 

In particular, the parametrical structure of asymptotic 
results for high K n can be easily understood in the coordi- 
nate representation since the characteristic radius of the 
potential is the Compton wave length of an electron h/m e c 
and the radius of atomic states is typically fin 2 /Zarac. 
Thus, the actual expansion is in n/K n , rather than just 
in 1/ K n - A similar situation is with low k„. Study of the 



by the Uehling potential |Zj 
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where the dispersion 'photon' mass 
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plays a role of the inverse Yukawa radius. Here and for 
other calculations in this paper we use relativistic units in 
which h = c = 1 , while for final results we restore c and h 
if necessary. 

The Lamb shift in the nonrelativistic approximation is 
of the form 

AE^(nl)= f drr 2 \R nl \ 2 Vu(r) 



(Za) 2 ^ F„i(k„) , (13) 
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where R n i(r) is the radial part of the Schrodinger wave 
function in a hydrogen-like atom 



<Pnim(r) = R„,i(r)Yi m (r/r) . 



(14) 



Applying the well-known analytic expression for R n i (r) to 
Eq. i(T3"|) and integrating over r, we obtain (see Eq. (f.9) in 
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After replacing the hypergeometric function by an explicit 
finite sum, we integrate over v and arrive at the following 
expression for F n i : 



The result for Fio 

F w (k) = -K 122 (k) + ^K 222 (k), (21) 

which follows from Eq. 1|16|) , is known in simpler terms and 
in particular in terms of elementary functions (see Eq.l[T)l). 
The general expression H21(l now presents a correction for 
any states in terms of elementary functions. Such an ex- 
pression is also very useful to derive various asymptotics 
once we find related asymptotics for F w (k). Another way 
of the F n i presentation as a single finite sum can be found 
in EH. 



3 Asymptotic behavior at large K n 

In the case of K n 3> 1 we can use asymptotics for the 
ground state function Fi (cf. [3] and Eq.lJBJ) 
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where the integrals 
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An expression for an arbitrary state can be also derived 
as an expansion over 1/ K n . Here we present few first terms 

(cf. uni) 



can be expressed in general in terms of the generalized 
hypergeometric functions (J3J 21 . Here we mainly follow 
our notation in [3], but the definition of the integral K 
(see also 0])) is different from the related integral I there. 
While in the nonrelativistic limit, when e = and the 
parameter c is integer, K a ^_ c (K) = I a ,b,c(K), in the rela- 
tivistic case with non-integer c the notation is K a ^_ c (K) = 

I a ,b,c+2e(K,e). 

We note that for integer a, b, c the result can be ex- 
pressed in terms of elementary functions. Using recursive 
relations (cf. [3|) 
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where ip(x) = r'(x)/ T(x) is the logarithmic derivative of 
the gamma function and 
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we express the correction for an arbitrary state through 
the expression for the ground state 



The results for the asymptotics of the difference (EjJ 
related to the Lamb shift are much simpler than the result 
for each level separately: 
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To test our calculations, we consider a limit ln« n 3> 1 
and find the leading logarithmic term within the effective 
charge approach with the help of a substitution 
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The result reads 
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The logarithmic contribution vanishes for the Lamb split- 
ting <P n ,i-i,i- The logarithmic results arc in agreement 
with the direct calculations above. 



We note that the expansion in (28) and (29) is effectively 
done in nj k„. Meanwhile, the leading term in (29) is sup- 
pressed by a factor of 1/n and the two first corrections are 
additionally suppressed by 1/n. 



4.2 Near circular states 



4 Asymptotic behavior at large n n and large n 

We note that the asymptotic coefficients depend on n and 
one may wonder about their behavior at high n. To study 
this we apply the well-known expansion for ip(z) at high 

z 

^ + f)=lnz + l-J-, + ... 
The result for the Uehling correction reads 
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where C = -ip(l) = 0.577215665... is Euler's constant. 
Certain simplifications are achieved once we do an as- 
sumption on a particular relation between values of I and 



4.1 Low-Z states 

An important feature of the result in Eq. (27) is that the 
parameter of expansion is rather n/n n than 1/k„. For in- 
stance, our explicit result for F n i at n > 1 and low I 
(I < n) is 
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We keep here the I dependence in the 1/n terms in order 
to derive a related result for the Lamb splitting 
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After studying n > 1 at low I, we turn to another case 
of n » 1 at low values of the radial quantum number 
n r = n — I — 1 ~ 1. In particular, n r = is related to the 
so-called circular state. In the limit of high K n and n we 
obtain 



"n,n— n r — 1 \K"n) 




In the same limit the specific difference related to the 
Lamb shift is 



^n,n— n, — 2,n— n, — 1 \^n) 
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The difference is suppressed by 1/n, as well as for low 
I, but, in contrast to Eq. l|^|l . there is no additional sup- 
pression. As a result, we see that the high-K„ expansions 
above (cf. 0111) are valid only in the case of k„ 3> n, 
which reduces the range of their applicability drastically. 
We consider the case of K n 3> 1, but not K n 3> n in Sect. El 



5 Asymptotics at low k 

In principle we are interested in high rather than in low 
k values, because the problem is related to muonic and 
exotic atoms. However, for high n, even for k>1 we can 
easily arrive at a situation when K n = n/n -C 1 and thus 
this region is of interest. 

The asymptotic behavior of F n i(n n ) at small values of 
K n was studied in j^j (see also (Hj). Various approaches can 
be used for that. One may start from our expression (|21|l 
(cf. JUKI) , taking into account that 
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or apply Eq. (|16|l with K a b c presented in terms of integral 
(|17fl or of generalized hypergeometric functions J5J. Actu- 
ally, the latter is the most straightforward way to obtain 
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a 1ow-k„ expansion. In case K n <C 1 the expansionreads 
(n + l)\K 2 j+ 2 j 1 (2Z + 4)!! 
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The first term of this expansion is obtained in [Hj and 
is in agreement with our expression. As one can see, the 
series is in fact over n ■ n n — k rather than n n . That sets 
a condition for applicability of the low-re„ asymptotics as 
K n <C 1/n. In particular, it means that the asymptotics 
Eq. l|3lfl) cannot be applied for Rydberg states even for the 
muonic hydrogen, i.e., for the smallest possible Z (Z = 1), 
where k ~ 1.5 and n n ~ 1.5 /n. 



6 High n asymptotic behavior 

We see that while we expand the generalized hypergeo- 
metric function in terms of either n n or 1/k„, the real pa- 
rameters of both expansions involve a factor of n directly. 
That is due to the increase of the coefficients of the K n - 
or 1/ n n - expansions with n which technically originates 
from the expansion of the factor 



K n VT 
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in the basic integral K a b c (K n ), while c = 2n. 

We note that a consideration of high n is not unreal- 
istic. For instance, in [neutral] antiprotonic helium for re- 
alistic levels HO] we find Z = 2, n ~ 30 > 1, k ~ 27 > 1, 
n n ~ 1. One of the reasons to study high-n states is that 
they very weakly interact with the nucleus, especially if a 
value of I is also high. Such an immunity to the nuclear- 
structure effects is an advantage from both theoretical and 
experimental point of view. Therefore and also because of 
simplifications in calculations we consider below circular 
or near circular states at n ^ 1. 



Once we know the general expression (0) in terms of 
3F2 (cf. |3]), we can consider in each order of the n n ex- 
pansion only terms leading in n (we did above a similar 
procedure to prove for the few first terms of series that 
the expansion is over nK n , and not over just n n ). 

Collecting the leading in n terms we arrive at the result 
in the limit n » 1, k„ « 1 and k<n 



^,2n(«n)^^B(n,| 



2F1 ( n,n; -; n 2 n 



2(nn n ) 2 F X ( n, n; -; K 2 n 



(36) 



We note that 
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at v <C n, where I v (z) is the modified Bessel function. For 
the latter one can apply the well-known explicit expres- 
sions for v — 1/2,3/2 and we arrive at the expression 
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To express the correction to energy F n i(K n ) in terms 
of the basic integrals Kb t 2 n (K n ) for near-circular states 
(n r = n — I — 1 ^C n) we need to transform the related 
coefficients in IJ16JI in the limit of high n. We note that 
the integral Kb,2n(K n ) does not depend on b in the lead- 
ing 1/rt-approximation and thus the I dependence of the 
correction comes from the I dependence of the coefficients 
of (|16|l . Eventually we find 
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To conclude this consideration we need to discuss the 
accuracy and validity of our derivation. It is valid for n n ^ 
1 and k ~ 1 and the corrections are of relative order 1/n. 
In the case of K <C 1 it is consistent with the leading term 
of the 1ow-k„ expansion l|33(l . 



6.2 The limit of high k„ for the near-circular states 



6.1 Limit of low k„ for the near-circular states 



For n n 3> 1 we also consider only near-circular states, for 

which b ~ n r = n — I — 1 <C n. We can rewrite Kf,,2n(i^n) 

The combination of the K abc integrals which actually en- in terms of the basic integrals as follows 

ters the equation for the vacuum-polarization energy shifts 
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We find that b -C n for the near-circular states, and, as 
long as we use l|Tfl|l . c = 2n for any state. 
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If b <C n, we can expand the exponential in the integrand 
and find 



Kb^ni^n) 



v 2 (l-v 2 /3) - 
dv — ttttt^- e 



(40) 
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that depends upon combination of parameters n n /n only. 
After a substitute of the variable t — l/Vl — v 2 in this 
integral we arrive at the result 
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Substituting the expression into the sum in Eo. l|16|) . 
and neglecting n — I as compared with n in coefficients of 
the sum, we obtain 
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Similarly to the previous subsection, we find that our 
derivation is appropriate for 1/k„ <C 1 and n/n n ~ 1 and 
the result has a relative uncertainty on the order of 1/n. 



6.3 Comparison of Iow-k and high-K asymptotics 

The region of highest interest is the one for high-K is 
when l/n n <C 1, but not n/n n <C 1, since for the op- 
posite situation (n/ K n 3> 1) we have already known the 
proper asymptotic form of the correction. We note that 
for the region of interest K n /n can be about unity or even 
larger (e.g., as in the case of l/n n ~ n^ 1 / 2 <C 1 and 
n/n n ~ n 1 / 2 ^> 1). In particular, if n n /n ^> 1 the result 
of the i-integration in Ea. l|42(l will mainly come from a 
narrow region (t — 1) <C 1. That means that we can im- 
prove the final result once we consider a complete series 
for the logarithm in Section RT2l and setting t = 1 for all 
terms except of the leading term of the expansion 



exp < 2n log ( 1 



1 + KnVl - V 2 



2n 



x e 



-irt*- 1 ) 



1 + K 
The final estimation reads for high n n 
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We can also rewrite the result for low K n (Eq. H38|) ') as 
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Fig. 1. Ratio of different F nt n-i asymptotics and its exact 
values for the circular state with n = 100: (a) - nK n <g 1 
(Eq.JSBl), (b) - ac„ < 1 (Eq.lEHJ), (c) - K n < 1 (Eq.gP), (d) 
- «n > 1 (Eq.ED), (e) - k„ » 1 (Eq.@3), (f) - «„/n > 1 
(Eq.J^OJ); the horizontal axis is related to F n , n -i- 



Comparing those two asymptotics for the circular-state 
correction, we find that the integral K2,2n(K n ) can be pre- 
sented as a product of a factor 



1 



2» 



(45) 



and a smooth function. The factor is varying in an ex- 
tremely broad region of K n : from low n n (k„ <C 1/n) to 
large K n (n n S> n), while the smooth function changes 
from being proportional to riT^I 2 to ln(«/n 2 ). Because of 
so smooth behavior we expect that the asymptotics with 
the explicit factor Eq. l|15|l can be successfully applied for 
a somewhat larger region, however, their accuracy there 
is unclear. Various asymptotics are compared to the exact 
result for n = 100 in Fig. ^ In particular we see that an 
explicit presentation of the factor of (k„/(1 + K n )) 2n re- 
ally improves agreement between the asymptotics and the 
exact solution. 



7 Other states 

Above we obtained the high-n asymptotic expressions in 
two specific regions of parameters where K n ~ 1/n and 
K n ~ n for the circular and the near-circular states only. 
We are also interested in finding asymptotics in these re- 
gions that are valid for low I. 

For low K n and I <C n we can use an approximate 
relation 

■^2(n-i-i),27i(«n) - ( 46 ) 

K 2 n n C dvv 2 (l-^pj (l-v 2 ) l+l cxp{-2n Kn VT^} , 
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Table 1. Asymptotics at x ^> 1 for the Uehling correction for 
the lowest s states. The correction is presented in terms of a 



dimensionless function F„i; AE(nl) 
F„i(k„) 



(o/tt) ((Za) 2 mc 2 /n 2 ) x 



F„o{x) 



-!ln(2z) + A±-£ 



3 



+ . 



: \ n (2x) + f - I + 



+ . 



■H2x) + f -§J + 



+ . 



: ln(2x) + ■ 



2vr i 24 36tt , 

~T + ~7I ZW T ■ 



■ ln(2x) + 



_ 5tt , 75 

90 2.t ' 2i^ 



and, neglecting I as compared with n in coefficients of the 
sum Eg. (|16|) . obtain an approximation 



F n l(Kn) 



(nn n ) 



2(1+1) 



dv v 2 1 - 



n-l-l 



exp < — 2nre„ \j 1 — v 2 > V^ 



i\(2l + i + l)\ 



l-v z 

1i 



(47) 



We can see that the asymptotic depends upon a combina- 
tion of parameters nre n , confirming the above-mentioned 
fact that it is the real parameter of expansion at low re n . 
In the other region corresponding to k„ ~ n we do 
not see a simple way to find a proper asymptotic form for 
low-! states. 



8 Summary 

Concluding, we have to briefly discuss corrections to the 
results derived. Since the parameter m/M (here M is the 
mass of the nucleus) in muonic, pionic and other exotic 
atoms is not as small as in conventional atoms, an impor- 
tant question is the accuracy of our results obtained in the 
external field approximation, i.e. in the limit m/M = 0. In 
antiprotonic atoms such effects are even more important 
than in the muonic case. The higher m/M corrections, 
which are quite important for muonic atoms, can be easily 
taken into account for the Lamb shift in the nonrelativis- 
tic approximation by substituting the mass of the orbiting 
particle m for the reduced mass ttlr = mM/(m + M) in 
Eq.JSJ) and re for kr = Zamn/m e . 

The functions F n i are presented above for arbitrary nl 
in a closed analytic form in various ways. Certain asymp- 
totics are also presented. The results at re n 3> 1 for the 
lowest states are summarized in Table ^ They are sim- 
ple and transparent. The Lamb shift result in Eq. (|16[) is 
obtained in the nonrelativistic approximation and is valid 
for any hydrogen-like atom as far as the relativistic cor- 
rections can be neglected. Some asymptotic results for the 
splitting of levels with Al = 1 at re n 3> n for some low ly- 
ing states are summarized in Table |2 

We studied the applicability of naive low-re and high- 
re asymptotics and found that the region where they are 



Table 2. Asymptotics of the Lamb-shift- induced differ- 
ence at i > 1 for the lowest states presented in terms 
of a dimensionless function $ n ab' AE(n,l) — AE(n,l') = 
(a /it) ((Za) 2 mc 2 /n 2 ) x <£„«>(«„) 



$nOi(aO 



4-7T 

3x y 



■1 + 



+ . 



2 , 1 
"15 + IP 

87T | 



•| + 



+ 



<Pni2(x) 



15 ~ 

47T , 

~pr + ■ 



-1 + 



+ . 



21 ~ 



+ . 



<£n23(a 



Stt , 
-■ZS + • 



■h + 



3 

7T 



$n34.{x) 



27 
40tt 



+ 



+ . 



10 10 10 





Fig. 2. Ratio of different F n>n -l asymptotics and its exact val- 
ues for the circular state with n — 30: (a) - nn n <C 1 f Eg, 1331 1. 
(b) - K„ < 1 (Eq.ldJO, (c) - K n < 1 (Eq.gD), (d) - Kn > 1 

(Eq.gTJ), (e) - re„ > 1 (Eq.^3), (f) - re„/n > 1 (Eq.^DJ); 
the horizontal axis is related to F nn -i. 



valid strongly depends on n. For high n we considered 

some additional asymptotics (see., e.g., Fig. [21 where the 

results are presented for a realistic value of n = 30). We 

found a sum of the leading terms for n»l expansions for 

both low-re and high-re cases. In particular, we found that 

most of the change by orders of magnitude of the Uehling 

correction in the circular states can be presented in terms 

/ \ 2 ™ . . 

of a simple factor I re„/(l + re„) ) which is multiplied by 

a smooth function. 

All results are obtained in the leading nonrelativistic 
approach and corrections due to that are of relative order 
(Za) 2 . The relativistic effects will be considered elsewhere. 
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